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(57) Abstract 

A method in a computer system for modeling a rhree-dimeasional structure of a model protein (112) is provided. In a preferred 
embodiment, rbe modeling is based upon a mree-dimeosioQal structure of a template protein (111) and an amino acid sequence alignment 
(1 13) of the model protein and the template protein. The proteins comprise a plurality of amino acids having backbone atoms and side chain 
atoms. For each amino acid in me model protein, when the template protein has an amino acid aligned with the amino acid of the model 
protein, the position of each backbone atom of the amino acid of the model protein is established based on the position of a topologically 
equivalent backbone atom in the aligned amino acid of (he template protein. The inter-atomic distance constraints for each pair of atoms 
with an established position is generated. Finally, the positioo of each atom in die model protein is set so that the inter-atomic distances 
are in accordance with the constraints. 
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Description 

METHOD AND SYSTEM FOR PROTEIN MODELING 

Cross-Ref erence to Prior Application 

This application is a continuation-in-part of U.S. 
Patent Application Serial No. 08/055, 0S0, filed April 28, 
1993, now pending. 

Technical Field 

This invention relates generally to a computer- 
implemented method and system for modeling the three- 
dimensional structure of a protein, and, more specifically, to 
protein modeling using the three-dimensional structure of a 
template protein and sequence alignment between the protein to 
be modeled and the template protein. 

Background Of The Invention 

The function of a protein is related to its three- 
dimensional structure. The three-dimensional structure of 
some proteins can be determined using X-ray crystallography 
or Nuclear Magnetic Resonance (NMR) . However, for other 
proteins, these methods cannot be used to determine the three- 
dimensional structure. When these methods cannot be used, 
computer-based protein modeling techniques have been used with 
some success. These protein modeling techniques use the known 
three-dimensional structure of a homologous protein to 
approximate the structure of another protein. 

In one such technique, the known three-dimensional 
structures of the proteins in a given family are superimposed 
to define the structurally conserved regions in that family. 
Among the members of a given family, there is considerable 
variation in the conformations of regions located between two 
consecutive structurally conserved regions, and thus, these 
regions are called the variable regions. These variable 
regions essentially contribute to the identity of a protein in 
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its family. However, the modeling of the variable regions has 
been unsatisfactory. 

Conventional homology modeling techniques have been 
used routinely to build models of proteases and antibodies. 
These techniques generally model the three-dimensional 
position of amino acids in structurally conserved regions by 
taking the Cartesian coordinates from homologous amino acids 
in a template protein with known three-dimensional structure. 
For the amino acids in the variable regions, these techniques 
take suitable loops from the Protein Data Bank (PDB) . (The 
PDB is a collection of protein information relating to the 
known structure of protein. The PDB is administered by 
Brookhaven National Laboratory.} Although these techniques 
are generally successful in modeling the structurally 
conserved regions for structurally undefined members of the 
family, these techniques have been unsuccessful in modeling of 
the variable regions. Since variable regions and structurally 
conserved regions of the model protein come from different 
protein structures, high energy short contacts are often found 
in the models. These high energy contacts are usually between 
inter-variable regions that are grafted from different known, 
protein structures. In cases where the sequence identity in 
the structurally conserved regions between the template and 
the model protein is weak, the interior amino acids are also 
susceptible to short contacts. Generally, these short 
contacts are removed by performing rotation around single 
bonds using interactive graphics, which is a tedious, and at 
times, an impractical procedure. Energy minimization has been 
used to relax strains in a model. However, the minimization 
procedure leads to structures that are trapped in local minima 
and relies entirely on the integrity of the starting 
structure . 

Certain proteins with very weak sequence identities 
fold into similar three-dimensional structures. For example, 
the three-dimensional structures of a number of helical 
cytokines fold in similar three-dimensional topology in spite 
of weak sequence homology. Members of helical cytokine family 
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not only show diversity in their disulfide topology but also 
the disulfide crosslinks are part of the variable regions. 
Besides the helical cytokines , there are other protein 
families with weak sequence identities and non-homologous 
disulfides in the variable regions. The prior homology 
modeling techniques produce unsatisfactory modeling of these 
families because of the absence of sequence homology in the 
structurally conserved regions and the presence of non- 
homologous disulfide crosslinks in the variable regions. 

Summary of the Invention 

It is an object of the present invention to provide 
a method and system for modeling the three-dimensional 
structure of members of a protein family/super family from the 
known ' three-dimensional structure of one or more members of 
that family. 

It is another object of the present invention to 
provide a method and system for modeling the three-dimensional 
structure of a protein by using position information of all 
topologically equivalent atoms in a template protein. 

It is another object of the present invention to 
provide a method and system for generating a sequence 
alignment between two proteins based on the resulting three- 
dimensional structure model. 

It is another object of the present invention to 
provide a method and system for generating a protein model 
that minimizes short contacts. 

It is another object of the present invention to 
provide a method and system for modeling proteins with weak 
sequence identity with a template protein. 

It is another object of the present invention to 
provide a method and system for modeling variable regions 
between two structurally conserved regions of a protein. 
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It is another object of the present invention to 
provide a method and system for modeling variable regions 
constrained by disulfide crosslinks of a protein. 

These and other objects, which will become apparent as 
the invention is more fully described below, are provided by a 
computer-implemented method and system for modeling a three- 
dimensional structure of a model protein. In a preferred 
embodiment, the modeling is based upon a three-dimensional 
structure of a template protein and an amino acid sequence 
alignment of the model protein and the template protein. The 
proteins comprise a plurality of amino acids having backbone 
atoms and side chain atoms. For each amino acid in the model 
protein, when the template protein has an amino acid align d 
with the amino acid of the model protein, the position of each 
backbone atom of the amino acid of the model protein is 
established based on the position of a topologically 
equivalent backbone atom in the aligned amino acid of the 
template protein. Then, the inter-atomic distance constraints 
for each pair of atoms with an established position is 
generated. Finally, the position of each atom in the model 
protein is set so that the inter-atomic distances are in 
accordance with the constraints. 

In another aspect, the present invention provides 
for a method and system for modeling a three-dimensional 
structure of a variable region of a model protein. The model 
protein has amino acids that have positions within a three- 
dimensional structure. The method receives relative 
positional information between pairs of amino acids. The 
relative positional information preferably includes the \y and 
<t> angle values between pairs of amino acids. The method then 
establishes a position for a first amino acid in the variable 
region. The method then for each amino acid pair in the 
variable region generates a model position for the amino acids 
based on the received relative positional information for the 
pair of amino acids. In a preferred embodiment the \|/ and $ 
angle values are collected for pairs of amino acid 
combinations in variable regions of a family of proteins. The 
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generated model position for the amino acids are preferably 
based on a random selection of the collected \y and $ angle 
values. 

Brief Description Of The Drawings 

Figure 1 is an overall block diagram of the input 
and output for the protein modeler. 

Figure 2 is an overall flow diagram of the protein 

modeler. 

Figure 3 is a flow diagram of a routine to generate 
a TEA position data structure for the model protein. 

Figure 4 is a routine to generate distance 
constraints for topologically equivalent atoms. 

Figure 5 is a routine to generate standard 
constraints for non- topologically atoms. 

Figure 6 is a flow diagram showing the DGEOM 

process. 

Figure 7 is a diagram illustrating the v and <fr 
angles between pairs of amino acids. 

Figure 8 is a diagram illustrating the amino acid 
angle matrix. 

Figure 9 is a flow diagram of a function that 
generates the amino acid angle matrix. 

Figure 10 is a flow diagram of a function to 
generate a variable region model for a model protein based on 
the amino acid angle matrix. 

Detailed Description Of The Invention 

The present invention provides, a computer-based 
method and system for modeling the three-dimensional structure 
of a protein (model protein) . In a preferred embodiment, the 
protein modeler generates the three-dimensional structure of 
the model protein based on a template protein with a known 
three-dimensional structure and an amino acid sequence 
alignment between the model protein and the template protein. 
Using the template protein and the sequence alignment, the 
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protein modeler generates a variety of inter-atom distance 
constraints for conserved regions and standard and chemical 
constraints for variable regions, and allows for the input of 
miscellaneous constraints (e.g., disulfide crosslinks). The 
protein modeler then generates a three-dimensional structure 
for the model protein using known techniques of distance 
geometry to ensure compliance with the constraints. 
Srinivasan and Sudarsanam, Protein Science 2:277-289 (1993), 
incorporated herein by reference, contains a further 
description of the present invention. 

The protein modeler generates the inter-atomic 
distance constraints for the model protein based on the 
position of topologically equivalent atoms (TEA) in the 
template protein. The protein modeler determines whether each 
atom in the model protein is topologically equivalent to an 
atom in the template protein. An atom of an amino acid of the 
model protein is topologically equivalent when the amino acid 
of the model protein aligns with an amino acid of the template 
protein and the amino acid of the template protein has a 
corresponding atom. If the aligned . amino acids are identical, 
then each atom of the model amino acid corresponds to the atom 
of the same name in the template amino acid. For example, the 
C- alpha atom of the model amino acid corresponds to the 
C-alpha atom of the template amino acid, and the C-beta atom 
of the model amino acid corresponds to the C-beta atom of the 
template amino acid. If the aligned amino acids are not 
identical, then each atom of the model amino acid may or may 
not have a corresponding atom in the template amino acid. 
Each backbone atom (N, C-alpha, C , and 0) has a corresponding 
atom. However, the correspondence between side chain atoms 
depends on the aligned amino acids. For example,, if the model 
amino acid is arginine and the template amino acid is alanine, 
then the C-beta atoms correspond. However, the C-gamma, 
C-delta, and other atoms of the arginine side chain have no 
corresponding atoms in the alanine side chain. 

The protein modeler sets the position of each 
topologically equivalent atom in the model protein to the same 
position as that of the topologically equivalent atom in the 
template protein. The protein modeler then calculates the 
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distance betw en each pair of atoms that have topologically 
equivalent atoms. These distances are the inter-atomic 
distance constraints for the topologically equivalent atoms in 
the model protein. In one embodiment/ these distance 
constraints can be relaxed by accounting for the overall 
homology between the model and template protein as discussed 
below. 

The protein modeler generates standard constraints 
based on the planarity of a side chain, the C-alpha and C-beta 
chirality, and C-alpha inter-atomic distances. 

In a preferred embodiment, the protein modeler 
assigns a helical conformation to each amino acid in the model 
protein and establishes a position for each atom. The protein 
modeler then uses well-known software packages to generate 
chemical constraints such as inter-atomic distances between 1- 
2 and 1-3 neighbors based on the connectivity of the model 
protein. The inter-atomic distance constraints are used to 
constrain the topology relating to non-topologically 
equivalent atoms (non-TEAs) . The selection of a helical 
conformation is arbitrary and used to establish position 
information for the standard software package. One skilled in 
the art would appreciate that one could generate such 1-2 and 
1-3 neighbor constraints without use of standard software 
packages, and, thus obviate the need to establish the position 
of each atom to generate these constraints. 

Once the various constraints are established, the 
well-known distance geometry software, DGEOM (discussed 
below) , is invoked to calculate the position of each atom in 
the model protein and thus generates the three-dimensional 
structure of the model protein. 

Figure 1 is an overall block diagram of the input 
and output for the protein modeler. The protein modeler 110 
is preferably implemented on a computer system. In a 
preferred embodiment, the protein modeler 110 inputs template 
protein structure 111, the model protein partial structure 
112, and the model/ template protein sequence alignment. Based 
on the input information, the protein modeler generates a 
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model protein structure 114. The format of the information 
defining the protein structures are in a format defined by the 
Protein Data Bank (PDB) . Table 1 illustrates the PDB protein 
structure format. 
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Table 1 contains an atom portion and a connectivity 
portion. The atom portion specifies for each atom in the 
protein a sequential atom number (iAtom), the atom name 
(Atom) , the containing amino acid name (AA) , the containing 
amino acid sequence number (iAA) , and the Cartesian 
coordinates (position) of the atom (X,Y,Z). The connectivity 
portion indicates the connectivity of the atoms in the 
protein. The connectivity specifies each of the (up to four) 
neighbors of an atom. 

Figure 2 is an overall flow diagram of the protein 
modeler. In step 201, the protein modeler generates an 
initial three-dimensional structure for the model protein. In 
a preferred embodiment, the protein modeler inputs model 
protein partial structure 112, which contains structure 
information in PDB format without position information. One 
skilled in the art would appreciate that the model protein 
partial structure 112 could be generated by the protein 
modeler based on the model protein amino acid sequence. The 
position information for each atom in the model protein is 
generated based on a helical conformation. This position 
information is used by DGEOM to determine chemical constraints 
between the atoms of the model protein. These chemical 
constraints are later used to constrain the position of non- 
topologically equivalent atoms. The arbitrary assignment of 
helical conformation to the peptide backbone and the choice of 
side chain rotomers does not impose any restriction on 
subsequent stages of the modeling procedure. In fact, non- 
helical conformations can also produce acceptable results. An 
acceptable helical conformation algorithm may be found 
Sundaram, K. and Srinivasan, S., Computer. Programs Biomed, 
"Computer Simulated Modeling of Biomolecular Systems," Vol. 
10, p. 29-34, 1979. As discussed above, this step is used to 
take advantage of predefined software packages that generate 
chemical constraints. However, the chemical constraints could 
be generated directly from the model protein amino acid 
sequence based on well-known chemical constraint principles 
such as 1-2 and 1-3 neighbor topology. 
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In step 202, the protein modeler generates a TEA 
position data structure for the model protein based on 
topologically equivalent atoms. The sequence alignment 
between the template and the model protein is used to generate 
the TEA position data structure for the protein model. The 
protein modeler retrieves the position information for each 
topologically equivalent atom in the template protein and 
stores the result in the TEA position data structure for the 
model protein. The protein modeler also indicates in the TEA 
position data structure which atoms have no topologically 
equivalent atom in the template protein. Table 2 illustrates 
a sample TEA position data structure (also known as a hybrid 
data structure because it contains information based on the 
model and template proteins). The TEA position data structure 
contains an entry for each atom in the model protein. Each 
entry contains the following fields. The derived field 
contains a f T f for a topologically equivalent atom and an 'M f 
for a non-topologically equivalent atom. The fields model 
amino acid number <miAA) , model amino acid name (mAA) , model 
atom number (miAtom) , and model atom name (mAtom) identify the 
atom of the model protein. The fields template amino acid 
number (tiAA), template amino acid name (tAA), template atom 
number (tiAtom) , and template atom name (tAtom) identify the 
topologically equivalent atom in the template protein. The 
fields X, Y, and 2 specify the position of a topologically 
equivalent atom. 
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In step 203, the protein modeler generates the 
distance constraints based on the TEA position data structure. 
The protein modeler determines the distance between each 
topologically equivalent atom. In step 204, the protein 
modeler generates standard constraints for non-topologically 
equivalent atoms. These standard constraints include 
planarity of side chains, 1-4 inter-C-alpha distances, and C- 
alpha and C-beta chirality. In step 205, the protein modeler 
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inputs miscellaneous constraints. These miscellaneous 

constraints allow for the methods of the present invention to 
be expanded to include other constraints as they become known. 
For example, a miscellaneous constraint could be distance 
constraints imposed. by disulfide crosslinks. In step 206, the 
protein modeler generates the three-dimensional structure for 
the "model protein based on the constraints and the helical 
three-dimensional structure of the model protein. The protein 
modeler preferably employs the well-known program DGEOM to 
generate the three-dimensional structure for the model 
protein. The DGEOM program is one of the various well-known 
programs that are based on the principles of distance 
geometry. Blaney et al., Quantum Chemistry Program Exchange. 
Program it 590 . (1989), Indiana University, Bloomington, 
Indiana, contains a description of DGEOM. A detailed 
explanation of the theory in practice of distance geometry can 
be found in Havel, T.F., Kuntz, I.D., and Crippen, G.M., Bull. 
Math Biol., "Theory and Practice of Distance Geometry," 
Vol. 45, pp. 666-720, 1983. 

Figure 3 is a flow diagram of a routine to generate 
a TEA position data structure for the model protein. This 
routine inputs the model and template sequence alignment and 
the three-dimensional structure of template protein, and 
outputs the TEA position data structure for the model protein. 
In step 301, the routine selects the next model and template 
amino acid from the sequence alignment. In step 302, if all 
the amino acids have already been selected, then the routine 
returns, else the routine continues at step 303. In step 303, 
if the model protein has an amino acid in the sequence 
position but the template protein has no amino acid, then the 
amino acid is part of a variable region and the routine 
continues at step 304, else the routine continues at step 3 05. 
In step 304, the routine adds an entry to the TEA position 
data structure for each atom of the model amino acid as a non- 
topologically equivalent atom and loops to step 301 to select 
the next amino acid in the sequence alignment. In step 305, 
if both the model and template proteins have an amino acid in 
the sequence alignment, then the routine continues at step 3 06 
to determine which atoms between the amino acids are 
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topologically equivalent, else the routine loops to step 301. 
In step 306, the routine selects the next atom in the model 
amino acid starting with a first atom. In step 307, if all 
the atoms in the model amino acid have already been selected, 
then the routine loops to step 301, else the routine continues 
at step 308. In step 308, if the template amino acid has a 
topologically equivalent atom to the selected atom of the 
model amino acid, then the routine continues at step 309, else 
the routine continues at step 310. In step 309, the routine 
adds an entry to the TEA position data structure for the 
selected atom as a topologically equivalent atom and loops to 
step 306. In step 310, the routine adds an entry to the TEA 
position data structure for the selected atom as a non- 
topologically equivalent atom and the routine loops to step 
306. 
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CODE TABLE 1 



void GenerateTEAPoaitionDataStructure ( ) 

( 

imAtom * 1; 
itAtom = 1; 

for i<«cAminoAcids; i++) 

{ tAA - Align (i] . tAA; 



mAA * Align [i] . mAA; 

if (tAA = '-' && mAA 
AddVarAAO ; 



if (tAA !» H -" «4 mAA !=• "- M > 
(for (j=l; j<-cmaxAtoms; 
if (PIMtmAAJ [j] — 1 
if (PIMftAA) tj) — 1) 

AddTEA ( imAtom++ , i tAtom++ ) ; 
else 

AddKTEA { iraAtom++ ) ; 

} 

J 



void AddVarAAO 

{ 

do 

AddNTEA ( imAtom++ ) ; 
while (Helical [imAtom] . iAA Helical [imAtom* 1] . iAA) 

I 

void AddNTEA (index x) 

{ 

Hybrid [x] .derived = "M"; 
Hybridtx] .imAA - Helical [x] .iAA; 
Hybrid [x] .mAA - Helical (x] . AA; 

Hybrid [x] .imAtom » Helical [x] .iAtom; 
Hybrid [x] .mAtom = Helical [x] .Atom; 



void AddTEA (index x, y) 

( 

Hybrid [x] .derived - "T H ; 

Hybridtx] .imAA - Helical(x) .iAA; 

Hybrid[x] .mAA * Helical[x] .AA; 

Hybridtx] .imAtom » Helical [x] .iAtom; 

Hybrid[x] .mAtom ■ Helical tx] .Atom; 

Hybridtx] .itAA - Templatty] .iAA; 

Hybridtx] .tAA » Templatty] .AA; 

Hybridtx] .itAtom « Templatty] .iAtom; 

Hybridtx] .tAtom = Templat [y] .Atom; 

Hybrid(x].x » Templatty] .x; 

Hybridtx). y ■ Templat [y] . y; 

Hybridtx]. z * Templat (y] . z; 



} 
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Code Table 1 contains "C" programming language 
pseudocode corresponding to the flow diagram of Figure 3. The 
data structure Align contains the protein sequence alignment. 
The data structure Helical contains data relating to the 
three-dimensional structure of a helical conformation of the 
model protein. However, the pseudocode does not use the 
helical position data. The data structure Template contains 
the data relating to the three-dimensional structure of the 
template protein. The data structure Hybrid contains the TEA 
position data structure. The fields names are the same as 
described above for the three-dimensional structure in PDB 
format and the TEA positional data structure. The pseudocode 
also uses the position identity matrix (PIM) defined in Table 
2. 

Table 2 is the position identity matrix for side 
chain atoms of amino acids. The rows of table 2 represent 
amino acids and the columns represent positional identity of 
the amino acids. Topologically equivalent atoms in different 
amino acids are aligned in columns with an- entry of 1. For 
example, if alanine is replaced by arginine in the model, the 
position of the arginine atoms — N, C-alpha, C, 0, and C-beta- 
-in the template protein are used. Although not shown, the 
position identity matrix also contains columns for each of the 
four backbone atoms. 
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n 

o 


(jl 








£.1 






Zl 


Z2 


Z J 


til 

HI 


Hi 


(51 


'SI 


Pi 


(8) 


(9) 


(101 


(lit 


(12) 


(13) 


(14) 


(IS) 


(16) 




Ala 1 


0 


0 


ft 


A 
U 


A 
U 


A 

u 


A 

u 


ft 

u 


A 
0 


0 


A 

0 


ft 

u 


Arg 1 


1 


0 


1 


A 


1 


A 


A 

0 


1 


A 

0 


0 


1 


1 


Asn 1 


1 


0 


1 


1 


A 

0 


0 


A 


A 


a 

0 


0 


A 

0 


A 


Asp 1 


] 


A 

0 


1 


1 


A. 


A 

u 


A 

0 


ft 


ft 


A 


ft 

0 


ft 

u 


Cys 1 




0 


0 


A 

0 


A 


A 


A 

u 


ft 

u 


ft 
(J 


ft 


ft 


ft 

u 


Glu 1 


1 


A 

0 


1 

1 


A 
U 


1 


1 


ft 

u 


A 

u 


ft 

u 


ft 

u 


ft, 


ft 


Gin 1 




U 


1 


ft 


1 

1 


1 
1 


ft 
u 


ft 
u 


ft 
u 


n 


ft 
w 


0 
\i 


Gly 0 




U 


n 
U 


n 
u 


A 
V 


n 
yj 


u 


n 
u 


ft 
u 


\j 


ft 

\J 


Q 


His 1 




0 


1 


1 


1 


1 


0 


0 


0 


0 


0 


0 


lie 1 




1 


1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


Leu 1 




0 


1 


1 


0 


0 


0 


0 


0 


0 


0 


0 


Lys 1 




0 


1 


0 


i 


0 


0 


1 


0 


0 


0 


0 


Met 1 




0 


1 


0 


1 


0 


0 


0 


0 


0 


0 


0 


Phe L 




0 


1 


1 


1 


1 


0 


0 


0 


0 


0 


0 


Pro 1 




0 


1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


Ser 1 




0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


Thr 1 




1 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


Trp 1 




0 


1 


1 


1 


1 


1 


0 


X 


1 


1 


0 


Tyr 1 




0 


1 


1 


1 


1 


0 


1 


0 


0 


1 


0 


Val 1 




1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 



Figure 4 is a flow diagram of the routine to 
generate distance constraints for topologically equivalent 
atoms. The routine inputs the TEA position data structure of 
the model protein and outputs distance constraints for the 
topologically equivalent atoms. In step 401, the routine 
selects the next topologically equivalent atom in the TEA 
position data structure, starting with the first topologically 
equivalent atom. In step 402, if all the topologically 
equivalent atoms have already been selected, then the routine 
returns, else the routine continues at step 403. In steps 403 
through 407, the routine loops determining the distance 
between the selected topologically equivalent atom and all 
other topologically equivalent atoms. In step 403, the 
routine chooses a next topologically equivalent atom after the 
selected topologically equivalent atom in the TEA position 
data structure. In step 404, if all the topologically 
equivalent atoms after the selected topologically equivalent 
atom have already been chosen, then the routine loops to step 
401 to select another topologically equivalent atom, else the 
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routine continues at step 405. In step 405, the routine 
calculates the distance between the chosen and selected 
topologically equivalent atoms. In step 406, the routine 
calculates the upper and lower bounds of the distance 
constraints. As discussed above, the upper and lower bounds 
can be relaxed based on the overall homology between the model 
and template proteins. In step 407, the routine outputs the 
distance constraints between the chosen and the selected 
topologically equivalent atoms and the routine loops to step 
403. Code Table 2 contains "C" programming language 
pseudocode for the flow diagram of Figure 4. 

CODE TABLE 2 



void GenerateTEAConstraints ( ) 

{ 

for (i*l; i<**nmAtoras ; i++> 

if {Hybridli] .derived «™ "T") 
for (j«i; j<*nmAtoms; 
if (Hybridtj] .derived - m T m ) 

(Calculate distance between atom i and j 
Output mAA and mAtom for i and j 

upper and lower bounds of distance 

); 

} . _ 

Figure 5 is a flow diagram of a routine to generate 
standard constraints for non-topologically equivalent atoms. 
The routine inputs the TEA position data structure and outputs 
the standard constraints for non-topologically equivalent 
atoms. In step 501, the routine selects the next non- 
topologically equivalent atom in the TEA position data 
structure, starting with the first non-topologically 
equivalent atom. In step 502, if all the ■ non-topologically 
equivalent atoms have already been selected, then the routine 
returns, else the routine continues at step 503. In step 503, 
if the selected non-topologically equivalent atom is part of a 
planar side chain, then the routine outputs a planar 
constraint in step 504 and continues at step 505. In step 
505, if the selected non-topologically equivalent atom is a C- 
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beta atom and the amino acid is isoleucine or threonine, then 
the routine outputs a C-beta chiral constraint in step 506 and 
continues at step 507. In step 507, if the selected non- 
topologicaily equivalent atom is a C-alpha atom and the amino 
acid is not glycine, then the routine outputs a C-alpha chiral 
constraint in step 508 and continues at step 509. In step 
509, if the selected non-topologically equivalent atom is a C- 
alpha atom, then the routine outputs the distance constraints 
for the next C-alpha atom in the peptide chain. The routine 
then loops to step 501 to select the next non-topologically 
equivalent atom. Code Table 3 contains "C" programming 
language pseudocode for the flow diagram of Figure 5. 

CODE TABLE 3 



void GenerateNTEAConstraints () 

[ 

for (i»l; KnmAtoms; i++) 

if (Hybrid [i] .derived - model) 

(if (Hybrid [i] .mAA is planar) 

Output ("plane", ids of atoms in plane); 

if (Hybrid. mAtom is CB £6 

Hybrid [i] .mAA is isoleucine or threonine) 
Output ( "chiral'*, ids of aide chain atoms, 
chiral volume) ; 

if (Hybridti] .raAtom is CA 4ft 

Hybridti] .mAA is not glycine) 
Output (" chiral id of CA atom, 

ids of side chain atoms, chiral volume) ; 

if (Hybridti] .mAtom is CA) 
Output (id of i atom, id of next CA atom, 
upper and lower distance bound) ; 

} 

} 
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Figure 6 is a flow diagram of the utilization of the 
DGEOM software package. In step 601, DGEOM initializes an 
inter-atomic distance matrix based on the helical three- 
dimensional structure of the model protein. DGEOM determines 
inter-atomic distances based on chemical constraints such as 
inter-atomic distances between 1-2 and 1-3 neighbors. 
Although DGEOM establishes these distance constraints for all 
atoms within the model protein, the distance constraints for 
topologically equivalent atoms are eventually overridden by 
distance constraints imposed by the template protein. In step 
602, DGEOM applies the distance constraints for the 
topologically equivalent atoms to the distance matrix. In 
step 603, DGEOM applies the standard constraints for non- 
topologically equivalent atoms to the distance matrix. In 
step 604, the routine applies miscellaneous constraints to the 
distance matrix. In step 605, DGEOM smoothes the upper and 
lower bounds in the distance matrix. In step 606, DGEOM 
generates the final three-dimensional structure of the model 
protein based on the smoothed distance matrix. 

The following example illustrates the methods of the 
present invention. Table 3 contains a sequence alignment for 
the murine interleukin-4 and human interleukin-4 proteins. In 
this example, the three-dimensional structure of the human 
interleukin-4 protein is known, and the three-dimensional 
structure of the murine interleukin-4 (the model protein) 
protein is to be modeled. In this example, a template protein 
is derived from the human interleukin-4 protein. The variable 
regions in the template protein are preferably expanded to 
allow for flexibility in modeling the transitions betw en 
conserved and variable regions. 

Continuing with the example, the protein modeler 
sets the position of each atom in the model protein with a 
topologically equivalent atom in the template protein to the 
position of the atom in the template protein. For example, 
the first aligned amino acids for these model and template 
proteins are both leucine. Since each atom in the model amino 
acid has a topologically equivalent atom in the template amino 
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acid/ the protein modeler sets the position of each atom. 
Once the position of each topologically equivalent atom in the 
model protein is established, then the protein modeler 
generates the inter-atomic distance constraints for each pair 
of atoms. The protein modeler then determines standard 
constraints for the non-topologicaliy equivalent atoms. For 
example, since the first amino acid, histidine, in the model 
protein does not align with an amino acid in the template 
protein, each atom of histidine is a non-topologically 
equivalent atom. Thus, a planar constraint, a C-alpha chiral 
constraint, and an inter-C-alpha distance constraint are 
generated. The protein modeler then inputs miscellaneous 
constraints that are generated by a user of the protein 
modeler. These constraints may typically be based on 
disulfide crosslinks or other empirically determined 
constraints. The protein modeler also generates chemical 
constraints for the non-topologically equivalent atom. 
Finally, the protein modeler applies all the constraints using 
the principles of the distance geometry to arrive at a three- 
dimensional structure of the murine interleukin-4 protein. 



TABLE 3 



MURIL4 : HGCDKNHLREIIGILNEVTGEGTPCTEMDVPNVLTATKNT 

HUMIL4 : -EAEAHKCDIT-LQEIIKTLNSLTEQKTLCTELTVTDIFAASKDT 

TEMPLT : LQEIIKTLNSLTEQKTLCTELTVTDIFAASKDT 

MURIL4 : TESELVCRASKVLRI FYLKHGK-TPCLKKNS SVLMELQ 

HUMIL4 : TEKETFCRAATVLRQFYSHHEKDTRCLGATAQQFHRHKQLIRFLK 

TEMPLT : TEKET FC RAATVLRQ FY S H HE K CL LIRFLK 

MURIL4 : RLFRAFRCLDSSISCTMNESKSTSLKDFLESLKSIMQMDYS - 

HUMIL4 : RLDRNLWGLAGLNSCPVKEADQSTLENFLERLKTIMREKYSKCSS 

TEMPLT : RLDRNLWGLAGLNSCPVKEADQSTLENFLERLKTIMREKYSKCSS 



In an alternate embodiment, the present invention 
provides a method and system for modeling a model protein 
that bears weak sequence identity to a template protein with a 
known three-dimensional structure. Proteins with weak 
sequence identity, but with similar three-dimensional 
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structure, are known as divergent proteins. In this 
embodiment, the template protein and the model protein are 
preferably divergent proteins. The positions of the backbone 
atoms of the structurally conserved regions of the model 
protein are allowed to deviate from the position of the 
topblogically equivalent atoms in the template protein. The 
deviation, called a "divergence factor, " is expressed by 
specifying that inter-atomic distances as ranges rather than a 
single value. These ranges allow the modeling to accommodate 
for differences between aligned amino acids that are not 
identical. That is, the position of the backbone atoms can be 
adjusted to accommodate the differences in the side chains. 

In one embodiment, all inter-atomic distances 
between pairs of atoms are allowed to deviate within the 
divergence factor. In particular, the inter-atomic distances 
between atoms in the same structurally conserved region are 
allowed to deviate {intra-region deviation), and the inter- 
atomic distances between atoms in different structurally 
conserved regions are allowed to deviate (inter-region 
deviation) . Alternatively, intra-region deviation is not 
allowed. By not allowing intra-region deviations, the 
backbone of each structurally conserved region is held rigid 
but is allowed to rotate relative to other structurally 
conserved regions. 

In this alternate embodiment, the protein modeler 
generates a three-dimensional structure of the model protein 
based on a template protein with a known three-dimensional 
structure, an amino acid sequence alignment between the model 
protein and the template protein, and an identification of the 
structurally conserved regions. Since the . sequence identity 
between the model protein and the template protein is assumed 
to be weak, accurate alignment and identification of 
structurally conserved regions can be achieved by the well- 
known techniques, such as, "secondary structure prediction" or 
"disulphide topology." 

To model a divergent protein, the protein modeler 
performs several modeling iteration with various divergence 
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factors. After each iteration, the number of distance 
violations for that iteration is recorded. A distance 
violation preferably occurs when a distance constraint cannot 
be satisfied to within a margin of error, which is .5 A. The 
margin of error can be increased or decreased to allow for 
flexibility in the modeling. In a preferred embodiment, the 
protein modeler starts with a divergence factor of 0 and 
increases the divergence factor by 1 A for each iteration. As 
the number of iterations increases, the number of violations 
typically tends to level off. The model generated when the 
number of violations levels off is typically selected as the . 
three-dimensional structure for the model protein. 

Modeling Variable Regions 

In another embodiment, the present invention 
provides a method and system for modeling the variable regions 
of a model protein using the \y and $ angles between amino acid 
pairs of a template protein or a family of template proteins. 
A variable region modeler generates a collection of y and $ 
angle values for each amino acid pair in a family of proteins. 
The variable region modeler creates a matrix of amino acid 
pairs, such that, each entry in the matrix is a list of all y 
and $ angle values that are in the variable regions of the 
family. The variable region modeler models the variable 
region by randomly selecting a \y and $ angle value from the 
matrix for each amino acid pair in the variable region and 
calculating the corresponding atomic position information 
(e.g. Cartesian coordinates) . 

Each variable region is bounded by a beginning 
structurally conserved region and an ending structurally 
conserved region. The variable region modeler generates the 
atomic position information for the amino acids in the 
beginning structurally conserved region from the v and 4> angle 
values in the corresponding structurally conserved region of a 
template protein. The variable region modeler then generates 
atomic position information for the amino acids in the 
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variable region. The variable region modeler successively 
selects each pair of amino acids in the variable region. For 
each pair of amino acids, the variable region modeler selects, 
preferably randomly, an \j/ and 4> angle value from the 
corresponding list in the matrix. The variable region modeler 
then calculates the atomic positions of the amino acids using 
the selected \y and 4> angle values. After processing the 
variable region, the variable region modeler generates the 
atomic positions for the amino acids in the ending 
structurally conserved region using the \j/ and <(> angle values 
of the corresponding structurally conserved region in the. 
template protein. Finally, the variable region modeler 
calculates the root mean square (rms) deviation between the 
ending structurally conserved region of the model protein and 
the template protein. The rms deviation indicates the degree 
to which the model for the variable region is consistent with 
the bounding structurally conserved regions. 

In a preferred embodiment, the variable region 
modeler iteratively generates models and rms deviations for a 
variable region. The variable region modeler selects the 
variable region model with smallest rms deviation. As the 
number of iterations is increased, the likelihood of 
generating successively smaller rms deviations also increases. 
In the preferred embodiment, the selection of \\? and 4> angle 
values is random, that is, a Monte Carlo method. 
Alternatively, other selection criterion can be used. For 
example, selection of the vj/ and 4> angle values for each pair 
on the amino acids that follow and proceed the pair in the 
model protein. Also, as discussed above, amino acid angle 
matrix is preferably generated from the variable regions in a 
family of proteins. Alternately, the matrix can be generat d 
from proteins not in the same family. Also, a matrix can be 
generated for each corresponding variable region in a family 
of protein. Each matrix contains the y and <fr angle values for 
the corresponding variable region in each protein in the 
family. One skilled in the art would appreciate that other 
techniques for generating the matrix may be used. 
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After the atomic positions for the variable regions 
in a model protein are generated, then the positions may be 
used by the protein modeler as constraints to improve the 
modeling process as described above. 

Figure 7 is a diagram illustrating the y and $ 
angles between pairs of amino acids. The backbone atoms of an 
amino acid are planar, that is they are located within one 
plane. The y and <fr angles measure the inter-amino acid 
rotation between the C-alpha and C atoms and the N and C-alpha 
atoms, respectively. These angles specify the thre - 
dimensional structure of the. protein. 

Figure 8 is a diagram illustrating the amino acid 
angle matrix. The rows and columns of the matrix 801 
correspond to each of the 20 standard amino acid. Each entry 
of the matrix is a list of vy and <j> angle values for the 
corresponding amino acid pairs. For example, the entry 802 
contains the angle values for the amino acid pair asparagine 
followed by arginine. 

Figure 9 is a flow diagram of • a function that 
generates the amino acid angle matrix. The function inputs 
the three-dimensional structure of a template protein and 
updates the matrix based on the v|/ and 4> angles of the amino 
acids in the variable regions of the template protein. In 
step 901, the function identifies the variable regions in the 
template protein. All regions that are not secondary 
structural regions (which are identified in the PDB) are 
variable regions. In step 903 through 908, the function loops 
processing each variable region and within each variable 
region each pair of amino acids. In step 903, the function 
selects the next variable region in the template protein, 
starting with the first variable region. In step 904, if all 
the variable regions have already been selected, then the 
function is complete, else the function continues at step 905. 
In step 905, the function selects the next amino acid pair in 
the selected variable region, starting with the first. In 
step 906, if all the amino acid pairs in the selected variable 
region have already been selected, then the function loops to 
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step 903 to select the next variable region, else the function 
continues at step 907. In step 907, the function calculates 
the vy and <J> angles based on the position information of the 
selected amino acids. In step 908, the function adds the 
calculated \j/ and $ angle values to the list of angles in the 
matrix at the entry corresponding to the pair of the selected 
amino acids. The function then loops to step 905 to select 
the next pair of amino acids. 

Figure 10 is a flow diagram of a function to 
generate a variable region model for a model protein based on 
the amino acid angle matrix. The function outputs an rms 
deviation value for the generated model. In a preferred 
embodiment, this function is invoked multiple times, and the 
model with the smallest rms deviation is selected as the model 
for the variable region. In step 1001, the function sets the 
position information for the beginning structurally conserved 
region. The position information is generated from the y and 
$ angle values of the corresponding structurally conserved 
region in the template protein. Alternatively, the 
coordinates can be obtained directly from the three- 
dimensional structure of the template protein. In steps 1002 
through 1005, the function loops selecting each pair of amino 
acids in the variable region, selecting y and $ angle values, 
and setting position information for the atoms of the selected 
amino acids. In step 1002, the function selects the next 
amino acid pair in the variable region, starting with the 
first. In step 1003, if all the amino acid pairs in the 
variable region have been selected, then the function 
continues at step 1006 to process the ending structurally 
conserved region, else the function continues at step 1004. 
In step 1004, the function selects an v and an <fr angle value 
from the matrix at the entry corresponding to the selected 
amino acid pair. In a preferred embodiment, the V and $ 
angles are selected randomly from the list. Alternatively, 
the v and $ angles can be selected based on the probability of 
occurrence of the angles in the template protein or family of 
template proteins. Also, the v and 4 angle values can be 
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selected based on a probability that certain v and 4> angle 
values proceed or follow other v and <t> angle values. In step 
1005, the function calculates position information for each 
backbone atom in the selected pair of amino acids based on the 
angles and the calculated position information of the adjacent 
amino acid. The function also calculates the position 
information for the side chain atoms and loops to step 1002 to 
select the next pair of amino acids. In step 1006, the 
function determines the position information for the amino 
acids in the ending structurally conserved regions. The 
position information is calculated based on the calculated 
position of the last amino acid in the variable region and the 
y and <(> angles for each pair of amino acids in the ending 
structurally conserved region of the template protein. In 
step 1007, the function calculates the rms deviation between 
the ending structurally conserved region of the model protein 
and the corresponding region in the template protein. 

In alternate embodiment, the three-dimensional 
structure of a protein can be modeled de novo from the 
collected y and $ angle values. A matrix of \|/ and <fr angle 
values can be generated for all proteins in the PDB or for a 
family of proteins. The positional information c?n then be 
generated for each pair of amino acids in the model protein by 
selecting (e.g. randomly) \j/ and <fr angle values from the 
matrix. The compactness of the resulting three-dimensional 
structure can then be computed. The compactness can include a 
measurement of the volume of the three-dimensional structure, 
a measurement of the maximum inter-atomic distance, or a 
measurement of the density. The three-dimensional structure 
that is most compact can then be selected as the foled form of 
the model protein. Also, the v and <fr angle values can be 
combined with other constraints as described above to model a 
protein. 

Although the present invention has been described in 
terms of a preferred embodiment, it is not intended that the 
invention be limited to this embodiment. Modifications within 
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in the art, the scope of tha present invention is defined by 
the claims which follow. 
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Claims 

1. A method in a computer system for modeling a 
three-dimensional structure of a model protein, the modeling 
based upon a three-dimensional structure of a template protein 
and' an amino acid sequence alignment of the model protein and 
the template protein, the proteins comprising a plurality of 
amino acids, each amino acid having backbone atoms and side 
chain atoms, each atom in a three-dimensional structure having 
a position, the method comprising the steps of: 

for each amino acid in the model protein, when the 
template protein has an amino acid aligned with the amino acid 
of the model protein, establishing the position of each 
backbone atom of the amino acid of the model protein based on 
the position of a topologically equivalent backbone atom in 
the aligned amino acid of the template protein; 

generating inter-atomic distance constraints for 
each pair of atoms with an established position; and 

setting the position of each atom in the model 
protein wherein the inter-atomic distances are in accordance 
with the constraints. 

2. The method of claim 1 including the step of 
generating standard constraints for amino acids of the model 
protein having an atom with no topological equivalent atom in 
the aligned amino acid of the template protein. 

3. The method of claim 2 wherein the step of 
generating the standard constraints includes ' the step of 
generating a constraint specifying that .a side chain is 
planar. 

4. The method of claim 2 wherein the step of 
generating the standard constraints includes the step of 
generating a constraint specifying the chirality of a C-beta 
atom of an isoleucine or threonine amino acid. 
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5. The method of claim 2 wherein the step of 
generating the standard constraints includes the step of 
generating a constraint specifying the C-alpha chirality of an 
amino acid. 

6. The method of claim 2 wherein the step of 
generating the standard constraints includes the step of 
generating a constraint specifying inter-atomic distances 
between 1-4 neighbor C-alpha atoms. 

7. The method of claim 1 wherein the step of 
setting the position of each atom further includes the step of 
ensuring that the inter-atomic distance between a side chain 
atom and another atom in the same amino acid are in accordance 
with chemical constraints. 

8. The method of claim 7 including the step of 
generating chemical constraints based on 1-2 neighbor inter- 
atomic distances. 

9. The method of claim 7 including the step of 
generating chemical constraints based on 1-3 neighbor inter- 
atomic distances. 

10. The method of claim 1 wherein the step of 
establishing the position further includes the step of 
determining whether an atom of the template protein is 
topologically equivalent to an atom of the model protein. 

11. The method of claim 10 wherein the step of 
determining further includes the step of using a position 
identity matrix, the position identity matrix specifying 
topological equivalence of atoms in different amino acids. 

12. The method of claim 1 wherein the step of 
generating the inter-atomic distance constraints further 
includes the step of allowing for deviation of topology 
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between the model protein and the template protein based on 
the overall homology of the proteins. 

13. The method of claim 1 including the step of 
inputting miscellaneous constraints. 

14. The method of claim 13 including the step of 
inputting constraints based on disulfide crosslinks. 

15. The method of claim 1 including performing the 
steps for each chain of the model protein when the model 
protein has multiple peptide chains. 

16. The method of claim 1 including the step of 
deriving the template protein from multiple proteins. 

17. The method of claim 1 including the step of 
expanding variable regions of the template protein to allow 
for flexibility in the modeling of the transition from a 
conserved region to a variable region. 

18. The method of claim 1 including the step of, for 
each amino acid in the model protein, when the template 
protein has an amino acid aligned with an amino acid of the 
model protein, establishing the position of a plurality of 
side chain atoms of the amino acid of the model protein based 
on the position of topologically equivalent side chain atoms 
in the aligned amino acid of the template protein. 

19. The method of claim 18 including the step of 
generating standard constraints for amino acids of the model 
protein having an atom with no topological equivalent atom in 
the aligned amino acid of the template protein. 
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20. The method of claim 19 wherein the step of 
generating the standard constraints includes the step of 
generating a constraint specifying that a side chain is 
planar. 

21. The method of claim 18 wherein the step of 
establishing the position further includes the seep of 
determining whether an atom of the template protein is 
topologically equivalent to an atom of the model protein. 

22. The method of claim 18 wherein the step of 
generating the inter-atomic distance constraints further 
includes the step of allowing for deviation of topology 
between the model protein and the template protein based on 
the overall homology of the proteins. 

23. The method of claim 18 including the step of 
inputting miscellaneous constraints. 

24. The method of claim 18 including the step of 
inputting constraints based on disulfide crosslinks. 

25. The method of claim 18 wherein the step of 
setting the position of each atom further includes the step of 
ensuring that the inter-atomic distance between a side chain 
atom and another atom in the same amino acid are in accordance 
with chemical constraints. 

26. The method of claim 1 including the step of, for 
each amino acid in the model protein, when the template 
protein has an amino acid aligned with an amino acid of the 
model protein, establishing the position of each of the side 
chain atoms of the amino acid of the model protein based on 
the position of a topologically equivalent side chain atom in 
the aligned amino acid of the template protein. 
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27. The method of claim 25 wherein the step of 
setting the position of each atom further includes the step of 
ensuring that the inter-atomic distance between a side chain 
atom and another atom in the same amino acid are in accordance 
with chemical constraints. 

28. The method of claim 26 including the step of 
generating standard constraints for amino acids of the model 
protein having an atom with no topological equivalent atom in 
the aligned amino acid of the template protein. 

29. The method of claim 28 wherein the step of 
generating the standard constraints includes the step of 
generating a constraint specifying that a side chain is 
planar. 

30. The method of claim 26 wherein the steps of 
establishing the position further includes the step of 
determining whether an atom of the template protein is 
topologically equivalent to an atom of the model protein. 

31. The method of claim 26 wherein the step of 
generating the inter-atomic distance constraints further 
includes the step of allowing for deviation of topology 
between the model protein and the template protein based on 
the overall homology of the proteins. 

32. The method of claim 2 6 including the step of 
inputting miscellaneous constraints. 

33. The method of claim 26 including the step of 
inputting constraints based on disulfide crosslinks. 

34. A computer- implemented method for identifying a 
template protein for a family of topologically related 
proteins based on the modeling of a three-dimensional 
structure of a model protein, the proteins comprising a 
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plurality of amino acids, each amino acid having backbone 
atoms and side chain atoms, the template protein having a 
known three-dimensional structure, each atom having a position 
in a three dimensional structure, the method comprising the 
step of: 

for each amino acid in a model protein, when the 
template protein has an amino acid aligned with an amino acid 
of the model protein, establishing the position of each 
backbone atom of the amino acid of the model protein based on 
the position of topologically equivalent backbone atoms in the 
aligned amino acid of the template protein; 

generating inter-atomic distance constraints for 
each pair of atoms with an established position; 

setting the position of each atom in the model 
protein wherein the inter-atomic distances are in accordance 
with the constraints to generate a three-dimensional structure 
for the model protein; and 

accessing the conformance of the generated three- 
dimensional structure with the rules of protein folding. 

35. The method of claim 34 including the step of, 
for each amino acid in the model protein, when the template 
protein has an amino acid aligned with an amino acid of the 
model protein, establishing the position of a plurality of 
side chain atoms of the amino acid of the model protein based 
on the position of topologically equivalent side chain atoms 
in the aligned amino acid of the template protein, 

36. The method of claim 34 including the step of, 
for each amino acid in the model protein, .when the template 
protein has an amino acid aligned with an amino acid of the 
model protein, establishing the position of each of the side 
chain atoms of the amino acid of the model protein based on 
the position of a topologically equivalent side chain atom in 
the aligned amino acid of the template protein. 
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37. A computer- implemented method for obtaining a 
structurally accurate sequence alignment between a model 
protein and a template protein, the proteins comprising a 
plurality of amino acids, each amino acid having backbone 
atoms and side chain atoms, the template protein having a 
known three-dimensional structure, each atom having a position 
in a three dimensional structure, the method comprising the 
steps of: 

for a plurality of sequence alignments of the model 
protein and the template protein performing the steps of: 

for each amino acid in the model protein, when 
the template protein has an amino acid aligned with the amino 
acid of the model protein, establishing the position of each 
backbone atom of the amino acid of the model protein based on 
the position of topologically equivalent backbone atom in the 
aligned amino acid of the template protein; 

generating inter-atomic distance constraints for 
each pair of atoms with an established position; 

setting the position of each atom in the model 
protein wherein the inter-atomic distances are in accordance 
with the constraints to generate three-dimensional structure 
for the model protein; and 

accessing the conformance of the generated 
three-dimensional structure with the rules of protein folding. 

38. The method of claim 37 including the step of, 
for each amino acid in the model protein, when the template 
protein has an amino acid aligned with an amino acid of the 
model protein, establishing the position of a plurality of 
side chain atoms of the amino acid of the model protein based 
on the position of topologically equivalent side chain atoms 
in the aligned amino acid of the template protein. 

39. The method of claim 37 including the step of, 
for each amino acid in the model protein, when the template 
protein has an amino acid aligned with an amino acid of the 
model protein, establishing the position of each of the side 
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chain atoms of the amino acid of the model protein based on 
the position of a topologically equivalent side chain atom in 
the aligned amino acid of the template protein. 

40. A method in a computer system for modeling a 
three-dimensional structure of a model protein, the model 
protein having amino acids, the amino acids having atoms, each 
atom having a position in the three-dimensional structure, the 
method comprising the steps of: 

generating standard constraints for each amino acid 
in the model protein; 

generating an inter-atomic distance matrix for the 
model protein in accordance with the constraints; and 

setting the position of each atom in the model 
protein based on the distance matrix. 

41. The method of claim 40 wherein the step of 
generating the standard constraints includes the step of 
generating a constraint specifying the C-alptia chirality of an 
amino acid. 

42. The method of claim 40 wherein the step of 
generating the standard constraints includes the step of 
generating a constraint specifying inter-atomic distances 
between 1-4 neighbor C-alpha atoms. 

43. The method of claim 40 wherein the step of 
setting the position of each atom further includes the step of 
ensuring that the inter-atomic distance between a side chain 
atom and another atom in the same amino acid .are in accordance 
with chemical constraints. 

44. A method in a computer system for modeling a 
three-dimensional structure of a model protein, the model 
protein having amino acids, the amino acids having atoms, each 
atom having a position in the three-dimensional structure, the 
method comprising the steps of: 
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inputting miscellaneous constraints for each amino 
acid in the model protein; 

generating an inter-atomic distance matrix for the 
model protein in accordance with the constraints; and 

setting the position of each atom in the model 
protein based on the distance matrix. 

45. The method of claim 44 wherein the miscellaneous 
constraints are based on disulfide crosslinks. 

46. A method in a computer system for modeling a 
three-dimensional structure of a variable region of a model 
protein, the model protein having amino acids, the amino acids 
having positions within a three-dimensional structure, the 
method comprising the steps of: 

receiving relative positional information between 

pairs of amino acids; 

establishing a position for a first amino acid of 

the variable region; and 

for each amino acid pair in the variable region, 
generating a model position for the amino acids based on the 
received relative positional information for. the pair of amino 
acids . 

47. The method of claim 46 wherein the relative 
positional information includes v and $ angle values between 
pairs of amino acids. 

48. The method of claim 47 including the step of, 
for each combination of pairs of amino acids in variable 
regions of a family of proteins, collecting the vy and ♦ angle 
values for each pair of amino acids and wherein the step of 
generating a model position bases the model position on one of 
the collected vy and <(> angle values. 

49. The method of claim 4 8 wherein, the step of 
generating a model position bases the model position on a 
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randomly selected one of the collected v and $ angle values. 

50. The method of claim 4 6 including the step of 
generating a model position for the amino acids of the 
adjacent structurally conserved regions based on relative 
position information between pairs of amino acids. 

51. The method of claim 50 wherein the step of 
generating a model position for the amino acids of the 
adjacent structurally conserved region is based on the vj/ and <t> 
angle values in a corresponding region of a template protein. 

52. The method of claim 50 including the step of 
comparing the generated model position for the amino acids of 
the adjacent structurally conserved region to positions in a 
corresponding region in a template protein to indicate 
effectiveness of the modeling. 

53. A method in a computer system for modeling a 
three-dimensional structure of a variable region of a model 
protein, the model protein having amino acids, the method 
comprising the step of establishing positional information for 
the amino acids in the variable region based on y and 4> angle 
values between pairs of amino acids in a template protein. 

54. The method of claim 53 wherein the model protein 
has a structurally conserved region that is adjacent to the 
variable region, and including the steps of: 

establishing positional information for the amino 
acids in the adjacent structurally conserved region based on 
the established positional information of the amino acids in 
the variable region and based on v and $ angle values in a 
corresponding structurally conserved region of the template 
protein; and 

comparing the established positional information of 
the amino acids of the adjacent structurally conserved region 
to positional information of the corresponding structurally 
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conserved region of the template protein to measure the 
effectiveness of the modeling. 

55. The method of claim 53 including the step of 
collecting vy and $ angle values for pairs of amino acids in 
the template protein and wherein the step of establishing 
bases the positional information on the collected m> and 4> 
angle values. 

56. The method of claim 55 wherein the step of 
establishing includes the step of randomly selecting a 
collected vy and <fr angle value. 

57. The method of claim 55 wherein the step of 
collecting collects v and $ angle values for pairs of amino 
acids in a family of template proteins. 

58. A method in a computer system for modeling a 
three-dimensional structure for a variable region of a model 
protein, the protein having amino acids, the variable region 
having a corresponding beginning structurally conserved region 
and a corresponding ending structurally conserved region, the 
method comprising the steps of: 

collecting \y and 4> angle values for pairs of amino 
acids in a family of template proteins; 

generating three-dimensional positional information 
for the amino acids in the beginning structurally conserved 
region; 

generating three-dimensional positional information 
for the amino acids in the variable region based on the 
collected y and $ angle values and based on the generated 
positional information for the beginning structurally 
conserved region; 

generating three-dimensional positional information 
for the amino acids in the ending structurally conserved 
region based on the generated positional information for the 
variable region and based on positional information for a 



WO 94/25860 



40 



PCT/US94/04822 



corresponding ending structurally conserved region in a 
template protein; and 

comparing the generated positional information for 
the amino acids in the ending structurally conserved regions 
to positional information for the amino acids in the 
cortesponding structurally conserved region in the template 
protein to indicate correctness of the model. 

59. The method of claim 58 including the step of 
randomly selecting collected v and <(» angle values when 
generating the positional information for the amino acids in 
the variable region. 

60. The method of claim 58 including the step of 
repeating the steps of generating and comparing, and including 
the step of selecting, as the model of the variable region, 
generated positional information for the amino acids in the 
variable region when the generated positional information of 
the amino acids in the ending structurally conserved region 
most closely compares to positional information in the 
corresponding structurally conserved region in the template 
protein. 

61. A method in a computer system for modeling a 
three-dimensional structure of a model protein, the model 
protein having amino acids, the method comprising the step of 
establishing positional information for the amino acids in the 
model protein based on v and $ angle values between pairs of 
amino acids in a template protein. 

62. The method of claim 61 including the step of 
collecting v and $ angle values for pairs of amino acids in a 
template protein and wherein the step of establishing bases 
the positional information on the collected y and $ angle 
values. 

63. The method of claim 1 wherein the step of 
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generating inter-atomic distance constraints for each pair of 
atoms includes the step of including a divergence factor into 
the distance constraints. 

64. The method of claim 1 wherein the step of 
generating inter-atomic distance constraints for each pair of 
atoms includes the step of including a divergence factor into 
the distance constraints for inter-region distance and not 
including a divergence factor into the distance constraints 
for intra-region distances. 
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